Contribution of final state interaction to the branching ratio of B — > J/ipD 
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To testify the validity of the perturbative QCD (pQCD) and investigate its application range, one 
should look for a suitable process to do the job. B — > J/ipD is a promising candidate. The linear 
momentum of the products is relatively small, so that there may exist a region where exchanged 
gluons are soft and the perturbative treatment may fail, so that the non-perturbative effect would 
be significant. We attribute such non-perturbative QCD effects into the long-distance final state 
interaction (FSI) which is estimated in this work. We find that the contribution from the FSI to 
the branching ratio is indeed sizable and may span a rather wide range of 10" 6 ~ 10~ 5 , and cover 
a region where the pQCD prediction has the same order. A more accurate measurement on its 
branching ratio may provide important information about the application region of pQCD and help 
to clarify the picture of the inelastic rescattering (i.e. FSI) which is generally believed to play an 
important role in B decays. 

PACS numbers: 13.25.Hw, 13.75.Lb 



I. INTRODUCTION 

It is well known that B physics may provide an ideal 
field for testing all the existing theoretical frameworks, 
methods and searching new physics beyond the standard 
model (SM). The reason is that because there exists a 
heavy flavor, some approximations, such as the 1/Mq 
expansion can be adopted, so that the results of perturba- 
tive calculations are more reliable. When the study gets 
deeper, defects of such theoretical frameworks have been 
unavoidably exposed and demand further improving the 
standing theories. The most fundamental problem is how 
to properly evaluate the hadronic matrix elements which 
are fully governed by non-perturbative QCD. Thanks to 
the factorization, one can separate the non-perturbative 
QCD effects from the perturbative parts and the later one 
is calculable in terms of the field theory order by order. 
Based on this picture, various theories, such as the naive 
QCD factorization, pQCD (perturbative QCD) and the 
soft-collinear-effective theory (SCET) etc. are invented 
to calculate the processes where B mesons are involved. 

The pQCD has been proved to be a successful approach 
in B physics, namely most of the results obtained in this 
approach are consistent with data of the Babar, Belle 
and CLEO experiments. In this approach, the infrared 
divergence is properly dealt with by taking into account 
the contribution of the transverse momentum of quarks 
kx- In this picture, the non-perturbative part is included 
in the wavefunctions of the initial B meson and the pro- 
duced hadrons. Obviously, as one factorizes the pertur- 
bative part out and calculates the quark-level transition 
amplitudes, he must assume that all the constituents 
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which participate in the reaction are not far from their 
mass shells and moreover, all the internal lines (no mat- 
ter quark line or gluon line) must be hard enough, so that 
the perturbative calculation can make sense. A natural 
question would be raised, whether the pQCD framework 
is complete, even though its validity is supposed to be 
respected. By the asymptotic freedom of QCD at higher 
energy region, the perturbative approach works perfectly 
well, however, if one or two internal lines can reach a low- 
energy region where they are not sufficiently hard, one 
can be convinced that at this region the pQCD approach 
fails, or cannot result in reasonable values. If the inter- 
nal lines are soft, one can conjecture that at this region 
the non-perturbative QCD would dominate and it could 
be attributed into the so-called FSI, or the re-scattering 
sub-processes. To identify the application range of pQCD 
and testify its validity, we need to look for such processes 
where some internal lines can be soft. 

The process B — ► J/tp + D just serves for the purpose. 
The direct weak transition of B — > J/ip+D occurs via an- 
nihilation between b and u or VF-exchange between b and 
d (usually, we just name both of them as "annihilation"), 
and a pair of cc can emerge from a gluon splitting. Since 
the tub ~ 5.3 GeV, mj/^ ~ 3.1 GeV, and mo ~ 1-87 
GeV, and the linear momentum of the products in the 
CM frame of B meson is as small as 0.9 GeV, thus there 
exists a region where the gluon- line is soft. Thus one 
needs to include the contribution from the long-distance 
effects in the theoretical calculations as well as the short- 
distance effects which are contributed by the hard gluon 
lines. The strategy is following. There have been some 
calculations on the decay width for B — ► J/tp + D in 
terms of pQCD approach which we suppose to be the 
contribution of the direct transition from B into the fi- 
nal state J/tp + D, and then in this work, we calculate 
the contribution from long-distance effects to the rates, 
and then we urge our experimental colleagues to carry 
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out an more accurate measurement to testify the validity 
of the whole theoretical framework. 

There have been some experimental attempts along 
the line. The CLEO collaboration once reported a slow 
J/tp bump in the inclusive spectrum of B —* J/tjj + X 
[l| , which later was confirmed by Belle [|J and Babar [l| . 
These experiments indicate that there exists an excess 
in the momentum spectrum of the J/tp recoiling mass at 
about 2 GeV. The branching ratio of the excess is 6 x 
10~ 4 . Along with these experiments, different theoretical 
explanations have been suggested in Ref. [4j, [f| [(J . 

These theoretical works were focusing on the calcu- 
lations of direct transitions. Accompanying these theo- 
retical hypotheses [1, d, @| , theorist and experimentalist 
indeed began to study the branching ratio of B — > J/tpD. 
Assuming the intrinsic charm cc inside the B meson, 
Chang and Hou suggested that the branching ratio of 
B — > J /tp should reach an order of magnitude of 10 -4 Q. 
In 2002, by the collinear factorization approach, Eilam, 
Ladisa and Yang once calculated the branching ratio of 
inclusive B — > J/tp, and obtained it to be 7.28 x 10~ 8 
Q. The Babar and Belle collaborations reported a neg- 
ative result for searching B — > J/tpD decay. The upper 
limits on the branching fractions is set as 1.3 x 10~ 5 and 
2.0 x 10~ 5 for B° — > J/tpD° respectively corresponding to 
the Babar and Belle experiments @, [1] , which show that 
the assumption of the intrinsic charm cc inside the B me- 
son should be excluded. Later Li, Lu and Qiao reexam- 
ined the B — > J/tp in the framework of the pQCD fcy fac- 
torization, and predicted B[B° -> J/tpD] = 3.45 x 10~ 6 

1- 

As discussed above, besides the theoretical calculations 
on the decay width in terms of pQCD, one needs to 
take the FSI more seriously. In this work, we are go- 
ing to evaluate this contribution in terms of the hadronic 

loops ESElEim. Namel y' 

we consider several sub- 
processes such as B — > D(*'tt(D(*> p) — > J/tpD. Here we 
suppose the transition hamiltonian can be written as a 
sum 

i 

where Hi corresponds to both the quark-level and 
hadron-level hamiltonians and then 

(J/^D\H\B) = (J/t(>D\H quark \B) 

+ {J/4>D\H had \n){n\H' quark \B) + ... 

where H quar k and H' quark are the hamiltonian at quark 
level, but contribute to different states (for example 
J/tpD or D^n, D^*'p etc.) and the intermediate states 
|n) are the corresponding states with appropriate quan- 
tum numbers. Indeed (J/tpD\Hh a d\n) is just the inelastic 
re-scattering amplitude which should be evaluated. 

The above formulation indicates that the two parts 
should interfere, but the relative phase between the two 
parts (or several parts) is hard to determine because the 



different amplitudes are caused by different hamiltoni- 
ans and there (so far) is no any symmetry to associate 
them yet. To estimate the order of magnitude of such 
long-distance effects, we simply suppose the interferences 
among different modes are constructive. The first matrix 
element (n\H' quark \B) where |n) can be Dtt etc., can be 
evaluated reliably in terms of pQCD, since there is a suf- 
ficient phase space. 

Based on the idea, we re-evaluate the branching ratio 
of B — > J/tpD by considering the contributions from the 
hadronic loop effect to B° — ► J/tpD°. 

This paper is organized as follow. We present the cal- 
culation of Hadronic loop contribution for B° — ► J/tpD° 
in Sec. [Til Then we present the formulation about the 
factorization of B -> D ( *h(p) in Sec. HUJ In Sec. |IVJ 
the numerical result is given. The last section is a short 
conclusion and discussion. 



II. HADRONIC LOOP CONTRIBUTION FOR 

B° J/tpD° 

The diagrams which determine the hadronic loop ef- 
fects on the rate of B° — ► J/tpD° decay are depicted in 
Fig. [lj which can be divided into two groups. The fist 
group includes Fig. [TJ(a)-(d) and exactly corresponds to 
the left diagram of Fig. [5] which is depicted by a process 
at the quark level. Definitely, there are quark lines flow- 
ing from initial state hadron to the final state ones and 
therefore are the OZI allowed. Another group including 
only Fig. [1] (e) corresponds to the the quark-level process 
and is shown at the right diagram of Fig. [2] and obviously 
is an OZI forbidden diagram. 

According to the OZI rule, the contribution of Fig. [1] 
(e) is much suppressed comparing with that from Fig. 
[TJ (a)-(d). This fact can be confirmed by comparing the 
coupling constants of D^'D^J/tp and p-nJ/tp- Qfrr.jN, is 
about three orders larger than that of gn(*)D(*l j/v> 03- 
Thus we can safely ignore the contribution from Fig. [1] 
(e). 

The effective Lagrangians at the hadron-hadron ver- 
tices in Ref. [l4[ are 

£*dd = ig^DD^ (d^DD - Dd»D) , (1) 

+D*»{il> v d vl Dl-d lt !l> v D* v )} ) (2) 
C^D'D = -g.^ D > D £^ aP d^ (d a D* p D + Dd a D* p ) , (3) 
C D , D7r = ig D * D7I {pi&KD - Dd^irD;) , (4) 
£>d* D'tt = -g D . D . K e> a ' afi d M DZnd a Di, (5) 
Cdd p - ig DDp (Dp>'d tl D-d ll Dp>>D), (6) 
£d*d p = ~g D *D P ^ aP {Dd»p v d a Dl + dpDtdaPpD) . 

(7) 
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changed at t-channel, the resultant amplitude is 
Abs {a) 



d 3 pi d 3 p 2 , , 44 , . 

-(27T) d (p B ~Pi -P2) 



2 J (2^) 3 2E, i2;r; : 2/^ 

x [~«(P3 ^ Pl)a]* 2 V^ff.D*.D*-(iP2) 

^[^mi,.]. (9) 



1-^ + ^ 



raj,, J q 2 — m 2 



D~ 



The amplitudes of the mode £>° — > 
D-{ Pl )p+(p 2 ,e p ) -> J/^(p3,e 3 )5 (p4) where 
and D* + are exchanged respectively, read as 

1 /" d 3 pi d 3 p 2 



(2TT) 4 S\p B ~ Pl -p 2 ) 

Ml 



2 J {2tt) 3 2E 1 (2tt) 3 2E 2 
xA(B° -» D-p+)i 2 g^ DD e 3li [i(p 3 - Pi)' 1 - ipf 
xi 2 V2g_DBp[(-«)p4^ - (-*)(P3 



u£ , P2P2 \ 



<t~2 r 2 2 



2/2 2 J W ' "".DJ 



(10) 



and 




FIG. 2: Left-hand and right-hand sides diagrams respectively 
correspond to the OZI allowed and OZI forbidden diagrams 
for the hadronic loop diagrams for decay of B° -» J/^D°. 



Here £0123 = +1, tt = r ■ n, p = t ■ p. The definitions of 
the charmed meson iso-doublets are 



D T = (D°,D-), D=(D°,D+), 

D* T = (d*°.d*-), d* = (d*°,d* + ) . (8) 



We will present the values of the necessary coupling con- 
stants in the section for the numerical computations. 

With the above preparation, we can write out the de- 
cay amplitude involving contributions from the diagrams 
of Fig. [1] in terms of the Cutkosky Cutting Rules. For 
the process B° -» D~ fa)^ (p 2 ) -» J/^<J>3,e 3 )D (p 4 ) 
where as shown in Fig. [1] the vector-meson D* + is ex- 



Abs^ 

1 f d 3 pi d 3 p 2 



(2tt) S (pb — P\ — P2) 
2 J {2tt) 3 2E 1 {2tt) 3 2E 2 V ' KP y F ' 

xA(B° D~ p + )(-i)g^D* D e^ aP {-i)p^ v 

x[-i(ps -pi)a](-i)V2g D *D P e^ TX 

x(-i)(pa - pi)*(i)P2r( - gzp + -^jr) 

x(-.9 X A + ^)^^^ 2 te 2 ,^,]. (11) 
V mi,, Jo 1 — m/U. 



D* 



D ' 



For Fig. □] (c), the amplitude of B° -> 

£*-(pi,e D -+h + (p 2 ) -> J/^{p 3 ,e 3 )D°(p A ) where Z?*+ is 
exchanged at t-channel, is 

Abs (c) 

1 f d 3 pi d 3 p 2 



(27t) 4 <5 4 (pb-Pi-P 2 ) 



2 7 (2tt) 3 2£i {2tt) 3 2E 2 
xA(B° -» L>*"7r+)(-?; 2 )^ D - jD *[j(P3 - 2pi) • e 3 £W 

X - *(2p3 -Pl)i/C3A* + «(Pl +P3) M £ 3 l /] 



D*D^{iP2a) 



:) 



TO B* ' <l m B* 



(12) 



The amplitudes of 5° — > D* (pi, £_d*+)p + (P2, e?) — > 
J/V>(p3, e^)D (p4) where Z? + and D* + are exchanged re- 
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spectively are 



Abs^ 

1 f d 3 pi d 3 p 2 



2 J {2,f2E 1 (2,f2E^ )Hi{PB - pi ~ P2) 
xA(B° -► D*- p+){-i) gi , D * D E XT ^(~i)p 3X e 3T 
x(i)Pia(i 2 )V2g DDp [(-i)p^ - (-»)(p 3 



)( 



9 2 



2r„2 ^ 



(13) 



and 



1 /" d 3 pi d 3 p 2 



-(2tt) 4 5 4 (pb - Pi — P2) 



2 J (2tt) 3 2E 1 (2tt) 3 2E 2 
xA(B° -> .D^p+jH 2 ).^^^ - 2p x ) ■ e 35AT 
x - i(2p 3 - pi) T e 3 \ + i(pi + p 3 ) \6 3t ] 
x(-i)V2g D , Dp e a ^ x (-i){p 3 - Pl ) a (i)p 2u} 



9i 



D' 



l D* ' 1 



(14) 



D' 



In the above expressions, T\q 2 , mf] denotes the form fac- 
tor (FF), which reflects the structure effect at the effec- 
tive interaction vertices. In this work, following Ref. (l3j 
we take the monopole form for FF as 



F[q 2 ,m 2 ^ 



A 2 -m 2 
A 2 - q 2 J' 



(15) 



where the phenomenological parameter A can be param- 
eterized as 

A = m 2 + aA QCD . (16) 

rrij stands as the mass of the exchanged meson at t- 
channcl which is depicted in Fig. [TJ 

The decay amplitude of B° — ► J / ijjD via the hadronic 
loop diagrams is 

M[B° L>(*)-7r+(p+) -> J/ipD°] 
= Abs^ + Aba<- b -V + Abs^ 

+Abs {c) + Abs^-V + Abs {d - 2) . (17) 



III. B -> D ( *\(p) DECAYS IN THE 
FACTORIZATION APPROACH 

In this section, let us turn to study B — > D^*'ir(p). 
We can reliably apply the factorization which has been 
proved to be valid to all orders of the strong coupling 
constant in the heavy quark limit [HI, [l6| , to calculate 



the amplitude. The essential non-perturbative quantities 
are light meson (71", p) decay constants and B — ► D(D*) 
transition form factors. 

The decay constants for pseudoscalar (P) and vector 
(V) mesons are defined as 

(P(q)\A^\0) = -ifpq^ 
(V(q,e)\V p \0) = fvmye;. (18) 

where vector and axial vector currents are = qi"f^q2 
and A^ = 917^7592 respectively and e is the polarization 
vector of V. 

The B — > H (H = D, D* ) transition form factors are 
conventionally parameterized as in [l7| 



U)\\] > \B) = F l (T)\P,-^-rq^ 



+ 



(9 2 ){ 
P-q 



(19) 



(D*(e)\V^\B) = V( f ] Wp e**PV, (20) 
(D*(e)\A^\B) ^ l \2m D ,A {q 2 ) e —^-q ll 



+ {m B +m D ,)A 1 {q 2 )(e* 



" q 2 ' % 



~A 2 (q 2 



e -q 



rriB + VdD' 



P, 



P-q 



-<h> 



,(21) 



where P = p B + P d m, q — Pb ~ PdM an d P ■ Q 



m B ~ m D(* 



The decay amplitudes for B — > D^Tv(p) are given as 
A(B° -> D-n+) 

= ^V; b V udai {iU(m% m 2 D )F BD (ml)} , (22) 

A{B° -» D-p+) 
Gf 

= 7I 



( ' F V; b V udai {2f p m p F BD {w 2 ){e*-p B )) , (23) 



A{B° -> D*-ir+) 

= ^J=v; b V udai {2Um D ,A BD '{ml){e* D , -p B )) 



(24) 



A(B° -> £>*>+) 



V2 



V* b V u dai< - ifpUip 



(e* D , ■ s* p ){m B 



+m D *)A BU (m 2 ) 



2A BD '(m 2 p ) 



-ie liVaP e l f£ D ,p p p' D 



a 2V BD '(m 2 ) 



(m B + m D *) 



(25) 
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IV. NUMERICAL RESULT 



The relevant input parameters which are employed in 
this work include: m B o = 5279.4 MeV, m D o = 1864.5 



' 5 GeV~ 
5.84, 



MeV, mj/^p = 3096.9 
m D ,± = 2010.0 MeV, m„ 
MeV, G F = 1.16637 x 10~ 
41.6 x lO" 3 0; qd* = 

9il>DD = 9ipD'D' = 7.71, g^D'D 

9dd p = 2.52, g D , Dp = 2.82 GeV" 1 
and f p = 216 MeV 0. 

The Wilson coefficient a\ has been calculated up to the 
next-to-leading order [lj| and we take the value 



MeV, m D ± = 1869.3 MeV, 
139.6 MeV, m p ± = 775.5 
V ud = 0.974, V cb = 
go-D'Tr = 9.08 GeV-\ 
= 8.64, gu-D'p = 
03;/* = 132 MeV 



ai = 1.05. 



(26) 




The momentum dependence of the transition form fac- 
tors in eqs. (|22|) - ([25l) possess the pole structures [19( as 



F(d 2 ) = 



F(0) 



1 - a( + b( 2 



(27) 



with £ = q 2 /m 2 B . F(Q), a and b are obtained by fitting 
data and their values are shown in Table |TJ With the 



F 


F(0) 


a 


b 


^0 


0.67 


0.65 


0.00 


F BD 


0.67 


1.25 


0.39 


A BD* 


0.64 


1.30 


0.31 


A? D * 


0.63 


0.65 


0.02 




0.61 


1.14 


0.52 


yBD* 


0.75 


1.29 


0.45 



TABLE I: The values of F(0), a and b in the form factors of 
B -> D w [3. 



values given in Table HI one obtains F^ D {ml) = 0.67, 
Ff D {m 2 p ) = 0.69, A^ D \ml) = 0.64, Af D '(m 2 p ) = 0.64, 
A% D *(m 2 p) = 0.63, V BD *(m 2 p ) = 0.77, which will be ap- 
plied to the later numerical calculation. 

In Fig. [3J we show the dependence of the branching ra- 
tio of B° — > J/tpD a on the phenomenological parameter 
a in eq. I16I which spans a range a — 1 ~ 3. Furthermore, 
Table [TT| presents the branching ratio of B° — * J/tpD 
with some typical values of a. 



a 


1.0 


1.5 


2.0 


2.5 


3.0 


B[B" -> J/4jD u ] 
{xlQ- 5 ) 


0.13 


0.54 


1.4 


3.0 


5.2 



TABLE II: The branching ratio of B° 
ing to several typical values of a. 



J/ipD° correspond- 



FIG. 3: The variation of the branching ratio of B° -> J/ipD° 
with q = 1 ~ 3. 



V. CONCLUSION AND DISCUSSION 

In this work, we calculate the contribution of the FSI, 
i.e. inelastic rescattering processes to the branching ra- 
tio of B° — > J/ipD° and find that it spans a relatively 
wider range of 10 -6 ~ 10 -5 which is obviously larger 
than the theoretically predicted value 10~ 8 JB| and com- 
parable with the pQCD prediction of 10 -6 [9J], moreover, 
it is also noted that as a = 3 is taken, it can be close to 
the experimental upper bounds. We hope that the Babar 
and Belle collaborations will further investigate the pro- 
cess and carry out more accurate measurements. Then 
not only an upper bounds will be given, but also a rather 
precise branching ratio can be obtained. 

The significance of this investigation to the present the- 
oretical frameworks is obvious as discussed in the section 
of introduction. Since in the process the linear momen- 
tum of the final products is not large, there can exist 
a region where the exchanged gluon is soft and appli- 
cation of pQCD might fail. This region should be fully 
governed by the non-perturbative QCD effects which are 
not involved in the conventional pQCD calculations, even 
though phenomenological wavefunctions of the hadrons 
can partly cover such effects. Thus we consider the FSI 
effects as an additional contribution to that of pQCD 
evaluation. However, on another aspect, one cannot in- 
deed determine the range where pQCD fails and this is 
exactly the goal of this work. 

As noticed, there exist some phenomenological param- 
eters in our calculations on FSI effects, such as a or A in 
eq. (|16[) and other uncertainties which are coming from 
the employed data, therefore, we can only trust the re- 
sults to the order of magnitude. However, the largeness 
of the contribution of the FSI should draw our atten- 
tion because it may change the whole scenario. In fact, 
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an accurate measurement can provide an ideal field for 
testing validity of pQCD. Since the FSI can result in a 
branching ratio as large as 1CP 6 ~ 1CP 5 , there is a region 
where the pQCD prediction and the contribution of the 
FSI have the same order, thus it is hard to clearly iden- 
tify individual contribution from both mechanism, unless 
accurate measured data are available. On other aspect, 
under the assumption that the present pQCD calculation 
is trustworthy to a certain accuracy, it is not hopeless to 
determine their fractions because the two factories in- 
deed have ability to carry out such precise, but difficult 
measurements . 

More concretely, if the future measurement confirmed 
a smaller branching ratio of about 1CP 6 , then one would 
make a careful analysis to distinguish between the two 
kinds of contributions. Furthermore, if the data are basi- 
cally consistent with the pQCD predicted value, it is in- 
dicated that pQCD works well, even though there might 
be a range where application of pQCD is dubious. In 
other words, for that case, the range where pQCD fails, 
does not make dominant contribution and the whole the- 
oretical framework should cover a much wider application 
range than was expected. Then we have to re-adjust the 
input parameters for calculating the FSI or set a more 
stringent constraint on them. It would be very helpful 
for gaining knowledge on the FSI which plays important 
roles in many decay and production processes. 

By contraries, if the future measurements on B° — > 



J/ipD° confirm that the branching ratio is obviously 
larger than 1CP 6 , the fact would indicate that the re- 
gion where pQCD fails, is important and should be re- 
considered. In that case, we can conclude that one must 
be careful as he applies the pQCD to evaluate physical 
processes with low energy scales. And the FSI may be 
a possible solution to the discrepancy. If it is true, a 
byproduct would be that one can further investigate de- 
tails about the methods for calculation on the FSI effects 
and determine the concerned parameters, since the " con- 
tamination" from the direct process which is evaluated 
in terms of pQCD is relatively small. 

As a conclusion, we would urge our experimental col- 
leagues to make a more accurate measurement on this 
process because its significance to our theory is obvious. 
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